Multiscale velocity correlations in turbulence 



R. Benzi 1 , L. Biferale 2 ' 4 and F. Toschi 3 ' 4 
1 AIPA, Via Po 14, 00100 Roma, Italy. 
2 Dipartimento di Fisica, Universita di Tor Vergata, Via della Ricerca Scientifica 1, 1-00133 Roma, Italy. 
3 Dipartimento di Fisica, Universita di Pisa, Piazza Torricelli 2, 1-56126 Pisa, Italy. 

4 INFM-Umtd di Tor Vergata. 
(February 5, 2008) 



OO 

On 
ON 

C 

a 



> 

o 
o 

o : 
00 . 

ON ■ 

"c : 

I ' 

o : 
^ . 

o ■ 



X 



Multiscale correlation functions in high Reynolds num- 
ber experimental turbulence and synthetic signals are inves- 
tigated. 

Fusion Rules predictions as they arise from multiplicative, al- 
most uncorrelated, random processes for the energy cascade 
are tested. 

Leading and sub-leading contribution, in both the inertial and 
viscous ranges, are well captured by assuming a simple multi- 
plicative random process for the energy transfer mechanisms. 

In stationary turbulent flows, a net flux of energy es- 
tablishes in the inertial range, i.e. from forced scales, L, 
down to the dissipative scale r^. Energy is transferred 
through a statistically scaling-invariant process, which is 
characterized by a strongly non-gaussian (intermittent) 
activity. Understanding the statistical properties of in- 
termittency is one of the most challenging open problem 
in three dimensional fully developed turbulence. 

Intermittency in the inertial range is usually ana- 
lyzed by means of the statistical properties of velocity 
differences, 5 r v{x) — v(x) — v(x + r). In particular, 
in the last twenty years [1], overwhelming experimen- 
tal and theoretical works focused on structure functions: 
Sp(r) — ((8 r v(x)) p ). A wide agreement exists on the fact 
that structure functions show a scaling behaviour in the 
limit of very high Reynolds numbers, i.e. in presence of a 
large separation between integral and dissipative scales, 
L/r d — > 00: 
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The velocity fluctuations are anomalous in the sense that 
exponents do not follows the celebrated dimensional 
prediction made by Kolmogorov, = p/3. In fact, 
C(p) are observed to be a nonlinear function of p, which is 
the most important signature of the intermittent transfer 
of fluctuations from large to small scales. 
In order to better characterize the transfer mechanism, 
it is natural to look also at correlations among velocity 
fluctuations at different scales. Multiscale correlations 
functions should play in turbulence the same role played 
by correlation functions in critical statistical phenomena. 
Recently, some theoretical work [2,3] and an exploratory 
experimental investigation [4] have been devoted to the 
behavior of multiscale velocity correlations: 
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with rd < r < R < L. When the smallest among the 
two scales r goes beyond the dissipative scales, r^, new 
properties of the correlation functions (2) may arise due 
to the non trivial physics of the dissipative cutoff. From 
now on, we will mostly concentrate on correlation func- 
tions with both r and R in the inertial range. Moreover, 
in order to simplify our discussion, we will confine our 
analysis for the case of longitudinal velocity differences. 

Stochastic cascade processes are simple and well known 
useful tools to describe the leading phenomenology of the 
intermittent energy transfer in the inertial range. Both 
anomalous scaling exponents and viscous effects [1,5] can 
be reproduced by choosing a suitable random process for 
the multiplier, W(r, R), which connects velocity fluctua- 
tions at two different scales, R > r. 

The main finding of this letter is that experimental mul- 
tiscale correlations (2) are in quantitative agreement, for 
any separation of scale r / R, with the prediction one ob- 
tains by using a pure uncorrelated multiplicative process 
for the energy cascade. 

The main idea turns around the hypothesis that small 
scale statistics is fully determined by a cascade process 
conditioned to some large scale configuration: 



S r v(x) = W(r, R) ■ S R v(x) 



(3) 



where, requiring homogeneity along the cascade process, 
the random function W should depend only on the ratio 
r/R. Structure functions are then described in terms 
of the W process: S p (r) = C p ([W (r/L)] p ), where C p = 
((Slv(x)) p ) if the stochastic multiplier may be considered 
almost uncorrelated with the large-scale velocity field. 
Pure power laws arise in the high Reynolds regime: in 
this limit we must have ([W(£)] p ) - (^) c(p) . In the 
same framework, it is straightforward to give the leading 
prediction for the multiscale correlation functions (2): 
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which becomes in the hypothesis of negligible correlations 
among multipliers: 
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This expression was for the first time proposed in [2] 
and considered to rigorously express the leading behavior 
of (2) when r/R — » as long as some weak hypothesis 
of scaling invariance and of universality of scaling ex- 
ponents in Navier-Stokes equations hold. Let us notice 
that, beside any rigorous claim, expression (5) is also 
the zevo-th order prediction starting from any multiplica- 
tive uncorrelated random cascade satisfying ([VF(-^)] P ) = 
S p (r)/S P (R). 

In this letter we want to address three main questions: 

(i) whether the prediction (5) gives the correct leading 
behavior in the limit of large separation of scales r/R <~ , 

(ii) if this is the case, what one can say about sub-leading 
behavior for separation r/R ~ 0(1), (iii) what happens 
to those observables for which the "multiplicative pre- 
diction" (5) is incorrect because of symmetry reasons. 
Indeed, let us notice that for correlation like : 



F hg (r,R) = ((6 r v) (5 R v) q ) 



(6) 



the multiplicative prediction gives: 
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Such a prediction is wrong because, if homogeneity can be 
assumed, Si (r) = for all scales r. In this case prediction 
(5) does not represent the leading contribution. 

In this letter we propose a systematic investigation of 
(2) in high Reynolds number experiments and synthetic 
signals. The main purpose consists in probing whether 
multiscale correlation functions may show new dynamical 
properties (if any) which are not taken into account by 
the standard simple multiplicative models for the energy 
transfer. 

Experimental data have been obtained in a wind tun- 
nel (Modanc) with Re\ = 2000. The integral scale is 
L ~ 20 m and the dissipative scale is — 0.31mm. 
Synthetic signals are built in terms of wavelet decom- 
position with coefficients defined by a pure uncorrelated 
random multiplicative process [6]. Such a signal should 
therefore show the strong fusion rules prediction (5) and 
it will turn out to be an useful tool for testing how much 
deviations from (5), observed in experiments or numer- 
ical simulations, are due to important dynamical effects 
or only to unavoidable geometrical corrections. 

First of all, let us notice that for any 1-dimensional 
string of number (such as the typical outcome of labo- 
ratory experiments in turbulence) the multiscale correla- 
tions (2) feel strong geometrical constraints. In particular 
we may always write down "Ward-Identities" (WI): 



S P (R - r) = ([{v(x + R) - v(x)) - (v(x + r) - v(x))] p ) (7) 
= b(k,p)(-) k F k , p _ k (r,R), (8) 

k=0, P 

where b(k,p) = p\/[k\{p — k)l]. 
For example, for p = 2 we have 

2F 1A (r,R) = S 2 (r) + S 2 (i?) - S 2 (i?-r) ~ 

+ ° ©]■*<«> o 

where the latter expression has been obtained by expand- 
ing S 2 (i? — r) in the limit r/R — > 0. 
For p = 3we have 

S 3 (R - r) = S 3 (R) - S 3 (r) + 3F 2 ,i(r, R) - 3iq, 2 (r, R) 

The "Ward-Identities" will turn out to be useful for 
understanding sub-leading predictions to the multiplica- 
tive cascade process. One may argue that in geometrical 
set-up different from the one specified in (2) the same 
kind of constraint will appear with eventually different 
weights among different terms. 

The main result presented in this letter is that all 
multiscale correlations functions are well reproduced in 
their leading term, — > 0, by a simple uncorrelated 
random cascade (5) and that their sub-leading contri- 
bution, ~ 0(1); are fully captured by the geometri- 
cal constrained previously discussed, namely the "Ward- 
Identities" . 

The recipe for calculating multiscale correlations will 
be the following: first, apply the multiplicative guess for 
the leading contribution and look for geometrical con- 
straints in order to find out sub-leading terms. Second, 
in all cases where the leading multiplicative contribution 
vanishes because of underlying symmetries, look directly 
for the geometrical constraints and find out what is the 
leading contribution applying the multiplicative random 
approximation to all, non-vanishing, terms in the WI. 

Let us check the strong fusion rules prediction (5) 
for moments with p > l,q > 1. In Figure 1 we 
have checked the large scale dependency by plotting 
F p , q (r, R) / S p+q (R) ■ S P (R) as a function of R at fixed, 
r, for different values of p, q. 

The expression (5) predicts the existence of a plateau (in- 
dependent of R) at all scales R where the leading multi- 
plicative description is correct. 

^From Figure 1 one can see that, in the limit of large 
separation R — > L at fixed r, F PtQ (r, R) / S p+q (R) ■ S P (R) 
shows a tendency toward a plateau. On the other hand, 
there are clear deviations for r/R <~ 0(1). Such devia- 
tions show a very slow decay as a function of the scale 
separation. 

In order to understand the physical meaning of the ob- 
served deviations to the fusion rules (5), we compare, in 
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Figure 1, the experimental data against the equivalent 
quantities measured by using the synthetic signal. We 
notice an almost perfect superposition of the two data 
sets, indicating that the deviations observed in real data 
can hardly be considered a " dynamical effect" . 

Using the WI plus our multiplicative receipt for p = 4 
we quickly read that the leading contribution to i<2,2 is 
£)( r C(2)) . 0(_rC( 4 )-C( 2 )) 5 while sub-leading terms scale as 
(9(/( 4 )), and as 0(r^) ■ 0(R^-^). 
This superposition of power laws is responsible for the 
slowly-decaying correlations in Fig. 1. The result so far 
obtained, i.e. that both the experimental data and the 
synthetic signal show the same quantitative behaviour, is 
a strong indication that multiscale correlation functions, 
at least for p > 1, q > 1, are in good agreement with the 
random multiplicative model for the energy transfer. 

For multiscale correlations where the direct applica- 
tion of the random-cascade prediction is useless, like 
Fi i9 (r, R), we use the WI plus the multiplicative predic- 
tion applied to all terms, except the Fi iQ . One obtains 
the expansion: 
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which coincides when q — 1 with the exact result (9) 
using £(3) = 1. 

In Fig. 2 we show the experimentally measured Fi t 2 
and the fit that we obtain by keeping only the first two 
terms of the expansion in (10). The fit has been per- 
formed by imposing the value for the scaling exponents 
£(2), £(3) measured on the structure functions, i.e. only 
the coefficients in front of the power laws have been fitted. 
As one can notice, the fit works perfectly in the inertial 
range. Let us remark that the correlation changes sign 
in the middle of the inertial range, which is a clear indi- 
cation that a single power-law fit (neglecting sub-leading 
terms) would completely miss the correct behaviour. 

Next we consider the WI for p — 3. Due to the fact 
that Ss(r) ~ r in the inertial range, one can easily show 
that the WI forces F12 ~ F21. Therefore we can safely 
state that also correlation functions of the form F Pt i feel 
non trivial dependency from the large scale R, at variance 
with prediction given in [3] using isotropic arguments. 

Let us summarize what is the framework we have found 
until now. 

Whenever the simple scaling ansatz based on the uncorre- 
cted multiplicative process is not prevented by symme- 
try arguments, the multi-scale correlations are in good 
asymptotic agreement with the fusion rules prediction 
even if strong corrections due to sub-leading terms are 
seen for small-scale separation r/R ~ 0(1). Subleading 



terms are strongly connected to the WI previously dis- 
cussed, i.e. to geometrical constraints. In the other cases 
(i.e. Fi iq (r, R)) the geometry fully determines both lead- 
ing and sub-leading scaling. 

All this findings, led us to the conclusions that multi- 
scale correlations functions measured in turbulence are 
fully consistent with a multiplicative, almost uncorre- 
cted, process. 

Nevertheless, the strong and slowly-decaying sub-leading 
corrections to the naive multiplicative fusion rules pre- 
dictions are particularly annoying for any attempts to 
attack analytically the equation of motion for structure 
functions; in that case, multiscale correlations at almost 
coinciding scales are certainly the dominant contribu- 
tions in the non-linear part of the equations [3]. Indeed, 
as shown in an analytical calculation for a dynamical toy 
model of random passive-scalar advection [7] , fusion rules 
are violated at small scale-separation and the violations 
are relevant for correctly evaluating the exact behavior 
of structure functions at all scales. 

When the smallest distance r is inside the viscous 
length, one can use the approach of multiplicative pro- 
cesses with multiscaling viscous cutoff [8]. Namely, for 
the correlation D\ iq (R) = ((d x v) 2 (S R v) q ) one obtains: 



D hq (R) ~ (S R v)" 



Sv(r d )^ 



(11) 



where r d is the dissipative scale. In the multifractal inter- 
pretation we say: S rd v = (r d /R) h ■ 5 R v with probability 
P h (r d ,R) = {r d /Rf- D{ - h \ Following [8] we have: 



5v(r d ) ■ r d ~ l^—j S R v ■ r d ~ v. 



(12) 



Inserting the last expression in the definition of Di >g (i?), 
we finally have: 



dfi(h) (5 R v) q+2 R- 2 



RSrv 



2(h-l)+3-P(h) 

~ S q+3 (R) 



vR 



(13) 



where we have used the fact that the multifractal pro- 
cess is such that v((d x v) 2 ) — ► 0(1) in the limit v — > 0. 
Expression (13) coincides with the prediction given in 
[3]. The above computation are easily generalized for 
any {(d x v)P(5 R v)i). 

Finally, let us remark that the standard multiplicative 
process may not be the end of the story, i.e. the dynam- 
ics may be more complex than what here summarized. 
For example, one cannot exclude that also sub-leading 
(with respect to the multiplicative ansatz) dynamical 
processes are acting in the energy transfer from large to 
small scales. This dynamical corrections must be either 
negligible with respect to the geometrical constraints or, 
at the best, of the same order. 
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A possible further investigation of such issue would be 
to perform a wavelet analysis of experimental turbulent 
data. ^From this analysis one may hope to minimize 
geometrical constraints focusing only on the dynamical 
transfer properties. 

Other possible candidates to investigate the above 
problem are shell models for turbulence, where geomet- 
rical constraints do not affect the energy cascade mecha- 
nism. Work in both directions is in progress. 
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FIGURE CAPTIONS 

FIGURE 1: 

Experimental and numerical F P:q (r, R)/ S p+q (R)- S p (r) at 
fixed r and changing the large scale R. Circles correspond 
to p = 2, q = 2, diamonds to p = 4, q = 2 for the exper- 
imental data. Squares correspond to p = 2, q = 2 and 
triangles to p = 4, q = 2 for the synthetic signal. Small 
scale r is fixed to r = 16 in units of the Kolmogorov 
scale. The data for p = 2, q = 4 have been shifted along 
the vertical axis for the sake of presentation. 
FIGURE 2: 

Experimental Fi^(r, R) at fixed r = 16 and at varying 
R. The integral scale L ~ 1 x 10 4 rd. Let us remark that 
the observed change of sign in the correlation implies the 
presence of at least two power laws. The continuos line 
is the fit in the region r < R < L obtained by using only 
the first two terms in (10). 
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